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One of the important perturbative ("hard") probes of hot and dense QCD matter is 
the medium-induced energy loss of energetic partons, so called "jet quenching", which is 
predicted to be very different in cold nuclear matter and in QGP, and leads to a number of 
phenomena which are already seen in the RHIC data on the qualitative level. The inclusion 
of jet quenching and other important collective effects, such as radial and elliptic flows, 
in the existing Monte-Carlo models of relativistic heavy ion collisions is discussed. Some 
issues on the corresponding physical observables at RHIC and LHC energies obtained with 
HYDJET++ model are presented. 
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I. INTRODUCTION 

One of the basic tasks of modern high en- 
ergy physics is the study of the fundamental 
theory of strong interaction (Quantum Chromo- 
dynamics, QCD) in new, unexplored extreme 
regimes of super-high densities and temperatures 
through the investigation of the properties of 
hot and dense multi-parton and multi-hadron 
systems produced in high-energy nuclear colli- 
sions 

0,1113. Indeed ' Q CD 

is not just a quantum 
field theory with an extremely rich dynamical 
content (such as asymptotic freedom, chiral sym- 
metry, non-trivial vacuum topology, strong CP 
violation problem, colour superconductivity), but 
perhaps the only sector of the Standard Model, 
where the basic features (as phase diagram, phase 
transitions, thermalisation of fundamental fields) 
may be the subject to scrutiny in the labora- 
tory. The experimental and phenomenological 
study of multi-particle production in ultrarela- 
tivistic heavy ion collisions is expected to pro- 
vide valuable information on the dynamical be- 
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haviour of QCD matter in the form of a quark- 
gluon plasma (QGP), as predicted by lattice cal- 
culations. 

Experimental data, obtained from the Rel- 
ativistic Heavy Ion Collider (RHIC) at maxi- 
mum beam energy in the center of mass sys- 
tem of colliding ions = 200 GeV per nu- 
cleon pair, supports the picture of formation of 
a strongly interacting hot QCD matter ("quark- 
gluon fluid") in the most central Au+Au (and 
likely Cu+Cu) collisions 0, 0, 0, 0| • This appears 
as significant modification of properties of multi- 
particle production in heavy ion collisions as com- 
pared with the corresponding proton-proton (or 
peripheral heavy ion) interactions. In particu- 
lar, one of the important perturbative ("hard") 
probes of QGP is the medium-induced energy 
loss of energetic partons, so called "jet quench- 
ing" [1, 0|, which is predicted to be very different 
in cold nuclear matter and in QGP, and leads 
to a number of phenomena which are already 
seen in the RHIC data on the qualitative level, 
such as suppression of high-p^ hadron produc- 
tion, modification of azimuthal back-to-back cor- 
relations, azimuthal anisotropy of hadron spec- 
tra at high pt, etc. [h], EH- A number of ap- 
proaches to the description of multiple scattering 
of hard partons in the dense QCD-matter have 



LP. Lokhtin, L.V. Malinina, S.V. Petrushanko, A.M. Snigirev, I. Arsene, K. Tywoniuk 



been developed [H, Ha, HI M, M, Ha G3, K 
taking into account the interference pattern for 
the emission of gluons with a finite formation 
time (QCD analog of the Landau-Pomeranchuk- 
Migdal effect in QED [H El) Besides the ra- 
diative loss (which is supposed to be dominant 
mechanism of medium-induced partonic energy 
loss), collisional loss due to elastic scatterings 
were_also quan tified in different theoretical mod- 

" nil m m, E3, m, m, m m hi . 



els [2 

On the other hand, one of the most spec- 
tacular features of low transverse momentum 
( "soft" ) hadroproduction at RHIC are strong col- 
lective flow effects: the radial flow (ordering of 
the mean transverse momentum of hadron species 
with corresponding mass) and the elliptic flow 
(mass-ordered azimuthal anisotropy of particle 
yields with respect to the reaction plane in non- 
central collisions). The development of such a 
strong flow is well described by the hydrodynamic 
models and requires short time scale and large 
pressure gradients, attributed to strongly inter- 
acting systems [33| . 

The heavy ion collision energy in Large 
Hadron Collider (LHC) at CERN will be a fac- 
tor of 30 larger then that in RHIC, thereby al- 
lows one to probe new frontiers of super-high tem- 
perature and (almost) net-baryon free QCD [34|, 
H, [H, H3] • The emphasis of the LHC heavy ion 
data analysis (at ^fs = 5.5 TeV per nucleon pair 
for lead beams) will be on the perturbative, or 
hard probes of the QGP (quarkonia, jets, pho- 
tons, high-pr hadrons) as well as on the global 
event properties, or soft probes (collective ra- 
dial and elliptic flow effects, hadron multiplicity, 
transverse energy densities and femtoscopic mo- 
mentum correlations) . It is expected that at LHC 
energies the role of hard and semi-hard particle 
production will be significant even for the bulk 
properties of created matter. 



II. JET QUENCHING AND FLOW IN 
MONTE-CARLO MODELS OF 
RELATIVISTIC HEAVY ION COLLISIONS 

Ongoing and future experimental studies of 
relativistic heavy ion collisions in a wide range 
of beam energies require the development of new 
Monte-Carlo (MC) event generators and improve- 
ment of existing ones. Especially for experiments 



which will be conducted at LHC, due to very high 
parton and hadron multiplicities, one needs fast 
(but realistic) MC tools for heavy ion event sim- 
ulation [H, [3^, H3] • The main advantage of MC 
technique for the simulation of high-multiplicity 
hadroproduction is that it allows a visual com- 
parison of theory and data, including if neces- 
sary the detailed detector acceptances, responses 
and resolutions. A realistic MC event generator 
should include a maximum possible number of 
observable physical effects which are important 
to determine the event topology: from the bulk 
properties of soft hadroproduction (domain of low 
transverse momenta px < lGeV/c) such as col- 
lective flows, to hard multi-parton production in 
hot and dense QCD-matter, which reveals itself 
in the spectra of high-py particles and hadronic 
jets. 

In most of the available MC heavy ion event 
generators, the simultaneous treatment of jet 
quenching and collective flow effects is absent. 
For example, the popular MC model HIJING 41] 
includes jet production and jet quenching, but it 
does not include flow effects. The partonic energy 
loss MC models PYQUENjH, JEWEL Q- 
PYTHIA 0, YaJEM lljfand MARTINI gg] 
combine a perturbative final state parton shower 
with QCD-medium effects and simulate the jet 
quenching in various approaches, but without 
considering the soft component of final hadronic 
state. The event generators FRITIOF [471 ] and 
LUCIAE [HI] include jet production (without 
jet quenching), while some collective nuclear ef- 
fects (such as string interactions and hadron 
rescatterings) are taken into account in LUCIAE. 
Another MC model THERMINATOR [H in- 
cludes detailed statistical description of "ther- 
mal" soft particle production and can repro- 
duce the main bulk features of hadron spectra at 
RHIC (in particular, describe simultaneously the 
momentum-space measurements and the freeze- 
out coordinate-space data), but it does not con- 
sider hard parton production processes. There is 
a number of microscopic transport hadron mod- 
els (UrQMD [5fJ, QGSM [Bj, AMPT [H], etc.), 
which attempt to analyze the soft particle pro- 
duction in a wide energy range, however they 
also do not include in-medium production of high- 
er multi-parton states. Another heavy ion event 
generator, ZPC [53f, has been created to simu- 
late parton cascade evolution in ultrarelativistic 
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heavy ion collisions. From a physical point of 
view, such an approach seems reasonable for very 
high beam energies (RHIC, LHC). However the 
full treatment of parton cascades may require sig- 
nificant amount of CPU run time (especially for 
LHC). 



III. HYDJET++ MODEL AND ITS 
APPLICATIONS FOR RHIC AND LHC 



HYDJET++ event generator [54| includes 
detailed treatment of soft hadroproduction as 
well as hard multi-parton production, and takes 
into account medium-induced parton rescatter- 
ing and energy loss. The heavy ion event in 
HYDJET+- 1- is the superposition of two inde- 
pendent components: the soft, hydro-type state 
and the hard state resulting from multi-parton 
fragmentation. Note that a conceptually simi- 
lar approximation has been developed in 55|, 56( ]. 
HYDJET+- 1- model is the development and con- 
tinuation of HYDJET event generator [42 j. and 
it contains the important additional features 
for the soft component: resonance decays and 
more detailed treatment of thermal and chemi- 
cal freeze-out hyper surfaces 57, [58[. The main 
program HYDJET+- 1- is written in the object- 
oriented C++ language under the ROOT envi- 
ronment [5^|. The hard part of HYDJET++ is 
identical to the hard part of Fortran-written HY- 
DJET and it is included in the generator struc- 
ture as a separate directory. When the generation 
of soft and hard components in each event at 
given impact parameter b is completed, the event 
record (information about coordinates and mo- 
menta of primordial particles, decay products of 
unstable particles and stable particles) is formed 
as the junction of these two independent event 
outputs. 

The details on physics model and simulation 
procedure of HYDJET++ can be found in the 
corresponding manual [HI]. The main features of 
HYDJET++ model are listed only very briefly 
below. 



A. Hard multi-jet production 

The model for the hard multi-parton part of 
HYDJET++ event is the same as that for HYD- 



JET event generator, and it based on PYQUEN 
partonic energy loss model [42]. The approach 
to the description of multiple scattering of hard 
partons in the dense QCD-matter (such as quark- 
gluon plasma) is based on the accumulative en- 
ergy loss via the gluon radiation being associ- 
ated with each parton scattering in the expanding 
quark-gluon fluid and includes the interference ef- 
fect (for the emission of gluons with a finite for- 
mation time) using the modified radiation spec- 
trum dE/dl function of decreasing temper- 
ature T. The model takes into account radiative 
and collisional energy loss of hard partons in lon- 
gitudinally expanding quark-gluon fluid, as well 
as realistic nuclear geometry. 

The Fortran routine for single hard nucleon- 
nucleon sub-collision PYQUEN was constructed 
as a modification of the jet event obtained 
with the generator of hadron-hadron interactions 
PYTHIA_6.4 The event-by-event simula- 

tion procedure in PYQUEN includes 1) gener- 
ation of initial parton spectra with PYTHIA and 
production vertexes at given impact parameter; 
2) rescattering-by-rescattering simulation of the 
parton path in a dense zone and its radiative and 
collisional energy loss; 3) final hadronization ac- 
cording to the Lund string model for hard par- 
tons and in-medium emitted gluons. Then the 
PYQUEN multi-jets generated according to the 
binomial distribution are included in the hard 
part of the event. The mean number of jets pro- 
duced in an AA event is the product of the num- 
ber of binary NN sub-collisions at a given impact 
parameter and the integral cross section of the 
hard process in NN collisions with the minimum 
transverse momentum transfer p™ m . In order to 
take into account the effect of nuclear shadowing 
on parton distribution functions, the impact pa- 
rameter dependent parameterization obtained in 
the framework of Glauber-Gribov theory 6l[ is 
used. 



B. Soft "thermal" hadron production 

The soft part of HYDJET++ event is the 
"thermal" hadronic state generated on the chem- 
ical and thermal freeze-out hyper surf aces ob- 
tained from the parameterization of relativistic 
hydrodynamics with preset freeze-out conditions 
(the adapted C++ code FAST MC [13, HI]). 
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Hadron multiplicities are calculated using the ef- 
fective thermal volume approximation and Pois- 
son multiplicity distribution around its mean 
value, which is supposed to be proportional to 
the number of participating nucleons at a given 
impact parameter of AA collision. The fast soft 
hadron simulation procedure includes 1) genera- 
tion of the 4-momentum of a hadron in the rest 
frame of a liquid element in accordance with the 
equilibrium distribution function; 2) generation 
of the spatial position of a liquid element and its 
local 4-velocity in accordance with phase space 
and the character of motion of the fluid; 3) the 
standard von Neumann rejection/acceptance pro- 
cedure to account for the difference between the 
true and generated probabilities; 4) boost of the 
hadron 4-momentum in the center mass frame of 
the event; 5) the two- and three-body decays of 
resonances with branching ratios taken from the 
SHARE particle decay table 62]. The high gener- 
ation speed in HYDJET+- 1- is achieved due to al- 
most 100% generation efficiency of the "soft" part 
because of the nearly uniform residual invariant 
weights which appear in the freeze-out momen- 
tum and coordinate simulation. 



Validation of HYDJET++ with 
experimental RHIC data 



It was demonstrated in [57j, 158J that FAST 
MC model can describe well the bulk properties 
of hadronic state created in Au+Au collisions 
at RHIC at yfs = 200A GeV (such as particle 
number ratios, low-p^ spectra, elliptic flow coef- 
ficients V2(pT,b), femtoscopic correlations in cen- 
tral collisions), while HYDJET model is capable 
of reproducing the main features of jet quench- 
ing pattern at RHIC (high-pr hadron spectra and 
the suppression of azimuthal back-to-back corre- 
lations) [12]. Since soft and hard hadronic states 
in HYDJET+- 1- are simulated independently, a 
good description of hadroproduction at RHIC in 
a wide kinematic range can be achieved, more- 
over a number of improvements in FAST MC and 
HYDJET have been done as compared to earlier 
versions. A number of input parameters of the 
model can be fixed from fittingthe RHIC data to 
various physical observables [541 ]. 

1. Ratio of hadron abundances. It is 

well known that the particle abundances in 



heavy ion collisions in a wide energy range 
can be reasonable well described within 
statistical models based on the assumption 
that the produced hadronic matter reaches 
thermal and chemical equilibrium. The 
thermodynamical potentials (Tb = 0.0285 
GeV, JT S = 0.007 GeV, JT Q = -0.001, the 
strangeness suppression factor j s = 1, and 
the chemical freeze-out temperature T ch = 
0.165 GeV have been fixed in [53] from fit- 
ting the RHIC data to various particle ra- 
tios near mid-rapidity in central Au+Au 
collisions at </s = 200 A GeV (n~ [tt + , 
P/tt-, K-/K+, K-/TT-, p/p, A/A, A/A, 
3/3, <f>/K~, A/p, 3-/7T-). 

2. ~Low-pT hadron spectra. Transverse 
momentum px and transverse mass mx 
hadron spectra (vr + , K + and p with mx < 
0.7 GeV/c 2 ) near mid-rapidity at different 
centralities of Au+Au collisions at y/s = 
200 A GeV were analyzed in [5^]. The 
slopes of these spectra allow the thermal 
freeze-out temperature T th = 0.1 GeV and 
the maximal transverse flow rapidity in 
central collisions p™ ax (6 = 0) = 1.1 to be 
fixed. 

3. Femtoscopic correlations. Because of 
the effects of quantum statistics and final 
state interactions, the momentum (HBT) 
correlation functions of two or more par- 
ticles at small relative momenta in their 
c.m.s. are sensitive to the space-time char- 
acteristics of the production process on 
the level of fm. The space-time param- 
eters of thermal freeze-out region in cen- 
tral Au+Au collisions at ^/s = 200,4 GeV 
have been fixed in [58|] by means of fitting 
the three-dimensional correlation functions 
measured for tt + tt + pairs and extracting 
the correlation radii -R s ide) -Rout and i?i on g : 
Tf {b = 0) = 8 fm/c, At/ (6 = 0) = 2 fm/c, 
R f (b = 0) = 10 fm. 

4. Pseudorapidity hadron spectra. The 

PHOBOS data on 77-spectra of charged 
hadrons [63] at different centralities of 
Au+Au collisions at -y/s = 200A GeV have 
been analyzed to fix the particle densities 
in the mid-rapidity region and the maxi- 
mum longitudinal flow rapidity r/ max = 3.3 
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(Fig. ED left). Since mean "soft" and "hard" 
hadron multiplicities depend on the cen- 
trality in different ways (they are roughly 
proportional to A r part (6) and Nb in (b) re- 
spectively), the relative contribution of soft 
and hard parts to the total event multi- 
plicity can be fixed through the central- 
ity dependence of dN/drj. The correspond- 
ing contributions from hydro- and jet-parts 
are determined by the input parameters 
fif th = 0.06 GeV and pf n = 3.4 GeV/c 
respectively. 

5. High-pr hadron spectra. High trans- 
verse momentum hadron spectra (pt > 
2 — 4 GeV/c) are sensitive to parton pro- 
duction and jet quenching effects. Thus fit- 
ting the measured high-p-p tail allows the 
extraction of PYQUEN energy loss model 
parameters. We assume the QGP forma- 
tion time To = 0.4 fm/c and the number 
of active quark flavours Nf = 2. Then 
the reasonable fit of STAR data on high-pj- 
spectra of charged pions at different cen- 
tralities of Au+Au collisions at a/s = 200^1 
GeV [64] is obtained with the initial QGP 
temperature To = 0.3 GeV (Fig. [[J right). 

6. Elliptic flow. The elliptic flow coeffi- 
cient v 2 (which is determined as the second- 
order Fourier coefficient in the hadron dis- 
tribution over the azimuthal angle <p rel- 
ative to the reaction plane angle ipR, so 
that V2 = (cos 2(ip — i/)r))) is an impor- 
tant signature of the physics dynamics at 
early stages of non-central heavy ion col- 
lisions. According to the typical hydrody- 
namic scenario, the values V2(pt) at low-pT 
(< 2 GeV/c) are determined mainly by the 
internal pressure gradients of an expand- 
ing fireball during the initial high density 
phase of the reaction (and it is sensitive to 
the momentum and azimuthal anisotropy 
parameters 5 and e in the frameworks of 
HYDJET++), while elliptic flow at high- 
Pt is generated in HYDJET++ (as well 
as in other jet quenching models) due to 
the partonic energy loss in an azimuthally 
asymmetric volume of QGP. The values 
of 6 and e were estimated in [54[ for dif- 
ferent centrality sets by fitting the mea- 



sured by the STAR Collaboration trans- 
verse momentum dependence of the ellip- 
tic flow coefficient v 2 of charged hadrons in 
Au+Au collisions at y/s = 200 A GeV 
Note that the choice of these parameters 
does not affect any azimuthally integrated 
physics observables (such as hadron multi- 
plicities, r]- and pT-spectra, etc.), but only 
their differential azimuthal dependences. 



D. Simulations with HYDJET++ at LHC 

The heavy ion collision energy at LHC will 
be a factor of 30 larger then that at RHIC, 
thereby allows one to probe new frontiers of 
super-high temperature and (almost) net-baryon 
free QCD. The emphasis of the LHC heavy ion 
data analysis (at \/s = 5.5 TeV per nucleon pair 
for lead beams) will be on the hard probes of the 
QGP (quarkonia, jets, photons, high-pr hadrons) 
as well as on the global event properties, or soft 
probes (collective radial and elliptic flow effects, 
hadron multiplicity, transverse energy densities 
and femtoscopic momentum correlations). It is 
expected that at LHC energies the role of hard 
and semi-hard particle production will be signifi- 
cant even for the bulk properties of created mat- 
ter. HYDJET++ seems to be an effective simu- 
lation tool to analyze the influence of in-medium 
jet fragmentation on various physical observables. 

Figure [5] shows the pseudorapidity distribu- 
tion of charged hadrons and transverse momen- 
tum distribution of pions obtained with HYD- 
JET++ default settings (in particular, p™ m = 7 
GeV/c) for 5% most central Pb+Pb events. The 
estimated contribution of hard component to the 
total event multiplicity is on the level ~ 55% 
here, what is much larger than as compared with 
RHIC (~ 15%). Of course, this number is very 
sensitive to the parameter p™ m — minimal pt 
of "non-thermalized" parton-parton hard scatter- 
ings. For example, increasing the value p2? m up to 
10 GeV/c results in decreasing this contribution 
down to ~ 25%. 

Some applications of HYDJET++ for high- 
Pt studies at the LHC were discussed during this 
Symposium [66, 67]. In particular, the spectac- 
ular predictions of HYDJET++ are possible re- 
ducing the femtoscopic correlation radii [66|, [6^| 
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FIG. 1: The pseudo-rapidity distribution of charged hadrons (left) and the transverse momentum distribution of 
pions (right) in Au+Au collisions at i/s = 200A GeV for different centrality sets. The points are RHIC data, 
histograms are the HYDJET++ calculations (solid - total, dotted - hydro part, dashed - jet part). 
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FIG. 2: The pseudorapidity distribution of charged hadrons (left) and the transverse momentum distribution 
of pions (right) in 5% most central Pb+Pb collisions at y/s = 5500^4 GeV (solid - total, dotted - hydro part, 
dashed - jet part), p^ in = 7 GeV/c. 



and elliptic flow 63, [69|] in heavy ion collisions 
as one moves from RHIC to LHC energies due to 
the significant contribution of (semi)hard compo- 
nent to the space-time structure of the hadron 



emission source. 



IV. SUMMARY 

One of the important perturbative ( "hard" ) 
probes of hot and dense QCD matter is the 
medium- induced energy loss of energetic partons, 
so called "jet quenching", which is predicted to 
be very different in cold nuclear matter and in 
quark-gluon medium, and leads to a number of 
phenomena which are already seen in the RHIC 



data on the qualitative level. Another spectacu- 
lar feature of hadroproduction at RHIC support- 
ing the pattern of hot QCD matter formation 
is strong collective flow effect (radial and ellip- 
tic flows). In order to quantify these collective 
effects in a wide energy range (up to the LHC), 
one needs the realistic Monte-Carlo models. Such 
event generators should take into account a max- 
imum possible number of observable physical ef- 
fects which are important to determine the event 
topology: from bulk properties of soft hadropro- 
duction (domain of low transverse momenta pt) 
to hard multi-parton production in hot QCD- 
matter, which reveals itself in spectra of high-pr 
particles and hadronic jets. 

Among other heavy ion event generators, 
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HYDJET++ model focuses on the detailed sim- 
ulation of jet quenching effect basing on the par- 
tonic energy loss model PYQUEN, and also re- 
producing the main features of nuclear collective 
dynamics by fast (but realistic) way. The fi- 
nal hadronic state in HYDJET++ represents the 
superposition of two independent components: 
hard multi-parton fragmentation and soft hydro- 
type part. The soft part of HYDJET++ is the 
"thermal" hadronic state generated on the chemi- 
cal and thermal freeze-out hypersurfaces obtained 
from the parameterization of relativistic hydrody- 
namics with preset freeze-out conditions. HYD- 
JET++ is capable of reproducing the bulk prop- 
erties of heavy ion collisions at RHIC (hadron 
spectra and ratios, radial and elliptic flow, fem- 
toscopic momentum correlations) , as well as high- 
Y>t hadron spectra. HYDJET+- 1- is an effective 
simulation tool to analyze the influence of in- 
medium jet fragmentation on various physical ob- 



servables at the LHC. 
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